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ABSTRACT 

We present the analysis of Suzaku and XMM-Newton observations of the broad-line radio 
galaxy (BLRG) 3C 111. Its high energy emission shows variability, a harder continuum with 
respect to the radio quiet AGN population, and weak reflection features. Suzaku found the 
source in a minimum flux level; a comparison with the XMM-Newton data implies an increase 
of a factor of 2.5 in the 0.5 - lOkeV flux, in the 6 months separating the two observations. 
The iron K complex is detected in both datasets, with rather low equivalent width(s). The 
intensity of the iron K complex does not respond to the change in continuum flux. An ultra- 
fast, high-ionization outflowing gas is clearly detected in the XIS data; the absorber is most 
likely unstable. Indeed, during the XMM-Newton observation, which was 6 months after, the 
absorber was not detected. No clear roll-over in the hard X-ray emission is detected, probably 
due to the emergence of the jet as a dominant component in the hard X-ray band, as suggested 
by the detection above ~ 100 keV with the GSO on-board Suzaku, although the present data 
do not allow us to firmly constrain the relative contribution of the different components. The 
fluxes observed by the y-ray satellites CGRO and Fermi would be compatible with the putative 
jet component if peaking at energies E ~ 100 Me V. In the X-ray band, the jet contribution to 
the continuum starts to be significant only above lOkeV. If the detection of the jet component 
in 3C 111 is confirmed, then its relative importance in the X-ray energy band could explain 
the different observed properties in the high-energy emission of BLRGs, which are otherwise 
similar in their other multiwavelength properties. Comparison between X-ray and y-ray data 
taken at different epochs suggests that the strong variability observed for 3C 1 1 1 is probably 
driven by a change in the primary continuum. 
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1 INTRODUCTION 

The unified model of Active Galactic Nuclei (AGN) predicts that 
the ability of an accretion disk-black hole system to produce pow- 
erful relativistic jets is the main reason of the observed differences 
between radio loud (RL) and radio quiet (RQ) sources. Geomet- 
rical effects, through the relative inclination of an axisymmetric 
system with respect to the line of sight, can exp lain the majority 
of the differences betwe en the various subclasses d Antonuccil 1 9931 : 
lUrrv & Padovanilll995h . Although this model does work well to 
first order, it leaves still open a number of problems. 

In particular, a key question still to be fully answered is 
why powerful relativistic jets are produced only in 10 - 20% 
of AGN, and if and how this is related to the structure 
of the accretion flow in the RL and RQ sources. In other 



words, an often debated issue is how accretion and ejecta 
in AGN are linked, and how these me chanisms work un- 
der dif ferent physical co n ditions ( e.g., IBla ndford & Znaiekl 
19771: IBlandford & Paynd [l982l: ISikora. Stawarz. & Lasotal 



20071; iGarofald 120091: iGarofalo. Evans. & Sambrunj |20iq 



' — 1 iii — i II i ■ ■ » 

Tchekhovskov. Narayan, & McKinnev 2010; for a parallel with 

X-ray binaries, see iFender. Gallo. & Russelll l2010h . A closely 
related question concerns the evolution of the jets, and its link 
with the observational subclasses of RL. It is commonly accepted 
that the relativistic jet originates in the innermost regions, and as 
the angle between the jet axis and the line of sight decreases, the 
importance of its emission increases, due to beaming effects. This 
simple picture does not account for the whole RL phenomenology, 
different global properties be ing associated to sources with differ- 
ent radio m orphology (e.g., lHardcastle. Evans. & Crostonl [20071: 
|Palvll2009aHb| ). 
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A key observational challenge in RL AGN is to disentangle 
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the jet contribution from the disk emission, which could allow us 
to better understand what are the main differences in the nuclear 
regions between RL and RQ sources. Previous X-ray observations 
of RQ and RL objects established that these s ources exhibit spec- 
tra w i th subtle but significant differences (e .g., Zdziars ki & Grandil 
l200lMSambruna. Eracle ous. & Mushotzkv 2002). The X-ray emis- 
sion of RQ objects appears to be described at first order by a pri- 
mary power-law continuum, with features resulting from the repro- 
cessing of this primary continuum from cold and warm gas, like 
the iron K line complex at 6 keV, the reflection component, and of- 
ten ionized absorption and emission features. The main observed 
difference between RL and RQ sources is that the features due to 
reprocessin g in RL AGN app ear to be weaker compared to the RQ 
population (Ballantvne 2007, and references therein). 

To account for these properties, several scenarios have been 
suggested. Highly ionized accretion disks, different from the stan- 
dard, cold disk ty pical of Seyferts, could be obtaine d if the accre- 
tion rates are high l lBallantvne. Ross. & Fabianl 20021) . A difference 
in the inner accretion disks of BLRGs and RQs could result in 
a small solid angle subtended in RLs by the disk to the primary 
X-ray source dEracleous. Sambruna. & Mushotzky||2000h . On the 
other hand, assuming the same disk structure in both populations, 
dilution by non-thermal jet emission cou ld weaken the reprocessing 
features dGrandi. Urrv. & Maraschll 2002b . 

It is however important to notice that in the last few years 
XMM-/Vewto/1 and Suzaku are changing this simplistic view: ob- 
servations of large samples of RQ AGN are showing a large spread 
in the X-ray properties of this class, with a wider range of contin- 
uum slopes and with several cases of sources having little or no 
reflection. It is thus emerging that the distinction between the two 
classes is not so sharp, but rather that RL objects populate one end 
of the distribution for Seyferts and quasars (OSOs; [Sambruna et al.l 
2009, and references therein). 

Strong signature s of low- and high- velo city outflows are rather 
common in RQ (e.g., iTurner & Milled 2 009). Warm absorbers (out- 
flowing with velocities v ~ 100 - 1000 km/s) are detected in the 
soft X-ray spectra of more than half of local Seyfert gala xies (e.g., 
ICrenshaw. Kraemer. & Georgel2003l : iBlustin et ali 2005 ) . Their lo- 
cation at ~ 1 - lOOpc suggests a possible association with the op- 
tical Broad Line Region (BLR) or Narrow Line Region. Ultra-fast 
outflows, showing velocities v > 100 00 km/s, have been fo und to 
be present in ~ 40% of local Seyferts (Tombesi et al. 2010a). Their 
location on sub-pc scales suggests a direct association with accre- 
tion disk winds/outflows. 

Regarding RL objects, only in the last few years have 
sensitive and broad-band observations started to detect disk 
winds/outflows in Broad Line Radio Galaxies (BLRGs). Chan- 
dra and XMM-Newton observations of 3C 382 revea led the pres- 
ence o f a warm absorber w ith outflow-like properties dReeves et al.l 
120091 ; iTorresi et al.l |2010|) . while 3C 445 shows signatures of 
soft X-ray photoionized gas in emission and absorption (see 



Sambruna. Reeves, & Braitol2007l:lReeves et al.l20ld:lBraito et all 



201 ll) . Our Suzaku observations provided evidence for ultra-fast 
disk outflows in 3/5 BLRGs (blueshift velocity v out ~ 0.04 - 
0.15 c, and mass out flow rates comparable to the accretion rates; 
iTombesi et alj2010bl) . 

Differences in terms of the strength of the reflection features, 
as well as the presence or lack of warm absorbers or high- velocity 
outflows, in sources showing different radio properties is one of the 
natural consequences of the so-called "ga p-paradigm" ( Garofalo 
2009: lGarofaloT Evans. & Sambrunal boiO). Its main ingredient is 
the relative orientation of the disk and the black hole spin; the abil- 



ity of the source of producing powerful jets is related to retrograde 
systems. In particular, this model accounts for the weaker reflection 
features observed in BLRGs, a natural consequence of the larger 
gap between the innermost stable circular orbit of the disk and the 
black hole in retrograde sources than in prograde objects. At the 
same time, the model provides a simple interpretation for the pres- 
ence or absence of signatures of disk winds, depending on the size 
of the gap region and the accretion efficiency. 

An additional complication in the study of the properties of 
the accretion flow in RL and RQ, is the presence in the former of 
emission due to the jet, that can mask the thermal emission (directly 
observed and/or reprocessed) from the disk. From the observational 
point of view, only through a wide energy coverage and good sen- 
sitivity at medium-hard X-rays, can we attempt to disentangle the 
contribution of the jet and the disk emission. The request of both 
simultaneous broadband coverage and high sensitivity in the iron K 
line region has now been achieved with Suzaku. We then started a 
program aimed to observe with this satellite the brightest (and best 
studied) BLRGs. Here we present the X-ray view of the BLRG 
3C 111 obtained thanks to our Suzaku data and an archival long 
exposure with XMM-Newton; Table[T]summarizes the main global 
properties of the source. 

The paper is organized as follows. The main properties of the 
source and previous X-ray results are presented in Section|2] Sec- 
tions [3] and [4] describe the X-ray observations and data analysis. 
Our results are discussed in Section|5] and finally in Section[6]we 
summarize our work. Throughout this paper, a concordance cos- 
mology with H = 71 km s~' Mpc" 1 , fi A = 0.73, and Sl m = 0.27 
dSpergel et alfeooj [2007) is adopted. The energy spectral index, 
a, is defined such that F v oc v' a . The photon index is T = a + 1. 



2 THE BLRG 3C 111 

3C 111 is a nearby (z = 0.0485), well-studi ed BLRG with broad 
optica l lines [FWHM(Ha) ~ 4800 km s' 1 ; lEracleous & Halperr] 
2003]. It exhibit s a Fanaroff-Riley II ra dio morphology, with a 
single-si ded jet {Linfield & Perlevl 1 19841) showing superluminal 
motion dVermeulen & Cohenl Il994h . From the measured proper 
motion [p. = (1.54 ± 0.2)" yr -1 ] and the maximum angular size of 
the radio lobes (6 = 273": lNilsson et al.lll993l) . lLewis et alj J2005>) 
inferred a range of 10° < i < 26° (21° < ( < 26° unless the source is 
a giant radio galaxy) for the inclination angle of the jet. A recent es- 
timate of the jet inclination of i ~ 19° has been derived with a Very 
Long Baseline Array (VL B A) monitoring (form the MOJAVE pro- 
gram; iKadler et al.l2008l) . This program revealed also the presence 
of a variety of components in the jet of 3C 111: a compact core, 
superluminal jet components, recollimation shocks, and regions of 
interaction between the jet and its surrounding medium. 

From the relation w ith the bulge luminosity, 

iMarchesini. Celotti. & Ferraresd d2004l) estimated a black hole 
mass for 3C 111 of M BH = 3.6 x 10' M Q . New estimates, derive d 
from measurements of the Ha width dChatteriee et al. I l201ll) . 
range from 1.5 x 10 8 M B to 2.4 x 10 8 M , more than a factor 10 
lower than the values de rived by Marchesi ni. Celotti. & Ferraresd 
(2004). As suggested bv lChatteriee et alj d201 ll) . this difference is 
probably due to the different extinction adopted in the two papers. 

Several high-energy observatories targeted 3C 111. Simulta- 
neous RXTE and XMM-Newton observations dLewis et alj l2005t) 
revealed a rather flat continuum (T ~ 1.63 - 1.75, depending on 
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Table 1. Main global properties of 3C 111. 



RA (J2000) 


Dec (J2000) 


Redshift 


WH,gal 


FWHM[Ha] 






R g 




logL bo i/L Edd 








[10 20 cnT 2 ] 


[kms -1 ] 


[deg.] 


[1O 8 M ] 


[10 14 cm] 


[erg s" 1 ] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


04hl8m21.3s 


+38d01m36s 


0.0485 


29.1 


4800 


19 


36.0 


10.7 


45.56 


-2.11 












10 - 26 


1.5-2.4 


0.44-0.71 


44.40 - 45.26 
















1.8 


0.53 


44.68 


-1.69 



(1) Right Ascension. (2) Declinatio n. (3) Cosmological redshift . (4) N eutral Galactic absorption column density JKalberla et al. 2005). (5) Full Width 
Half Maximum of the optical line jEracleous & Halpern 1994, 2003). (6) Inclinati on angle of the jet, from radio observations iKadler et alj 2008; 
Eracleous, Sambruna, & Mushotzky 2000); in the second row, the range inferred bv lLewisetal] j2005h is reported. (7) Black-h ole mass: in the first 
row, from Marchesini, Celotti, & Ferrarese 12004); in the second row, the range recently proposed by IChatteriee et afl 4201 ll) . which intermediate 
value is reported in the third row. (8) Schw arzschild radius; in each row, the valu e inferred from the corresponding black hole mass. (9) Bolometric 
luminosity: i n the first row, e s timate d by iMarchesini. Celotti. & Ferrarese 12004]) from the optical luminosity; in the second row, the range recently 
proposed by Chatte riee et alj 4201 lh . which intermediate value is reported in the third row. (10) Eddington ratio, from the black hole mass and 
bolometric luminosity reported in columns (7) and (9), respectively. 



the adopted model) and an extremely weak reflection component], 
R < 0.3. To explain the broad residuals found in the iron K energy 
range, the authors assumed either reprocessed emission (i.e., reflec- 
tion component and broad iron K emission line) from a truncated 
accretion disk, or a partial covering absorber with high column 
density (7Y H ~ 10 23 cm 4 ); although the two parametrizations were 
not distinguishable from the fitting point of view, the latter model, 
though less complex, was disfavoured by the authors on physical 
grounds (mainly in view of the small inclination angle of the jet). 
Regardless of the continuum model used, a narrow Fe Ka line aris- 
ing in a distant reprocessor is required, with equivalent width (EW) 
~ 20-30eV. 

A BeppoSAX observation of 3C 111 confirmed that a flat con- 
tinuum is present (r ~ 1.58), with only upper limits for the re- 
flection component and the Fe K a line (EW < 72 eV, R < 0.3; 
iGrandi. Malaguti. & Fioccm1l2006l ) . However, a r e-analysis of the 
same BeppoSAX data presented in lDadinj J2007h found a steeper 
continuum (T ~ 1.75) and a higher upper limit for the reflection pa- 
rameter, R < 2.25; the BeppoSAX data provides only a lower limit 
to the high-energy cutoff, E c > 82 keV. 

The 6-years of m onitoring with RXTE recently presented in 
IChatteriee et alj J20 1 lh highlights the long time-scale variability of 
3C 111. The flux shows a range of variability F2-iokcV ~ 2xl0 -11 - 
8 x 10~" ergs ctrT 2 s -1 (i.e., a factor of ~ 4), clearly correlated 
with the variation in the X-ray photon index. From their analysis 
of the variations in the X-ray continuum and Fe Ka emission line, 
the authors conclude that the iron line is generated within 90 light- 
days of the source of the X-ray continuum. The latter can be either 
the corona or the base of the jet; in both cases, the RXTE data are 
consistent with the standard paradigm of X-ray emission dominated 
by reprocessing of thermal photons produced in the accretion disk; 
no strong contribution from the jet is observed. 

A lOks Chandra exposure detected an excess of X-ray emis- 
sion in 3 radio knots out of the 4 present al ong the single sid ed jet 
structure, as well as in its terminal point dHogan et al.ll201 lh . The 
emission is ascribed by the authors to inverse Compton scattering 
off of cosmic microwave background photons (IC/CMB). Thus, de- 
pending on the adopted model for the jet structure from pc to kpc 
scales (e.g., jet bending and/or deceleration, or neither of them), the 



1 The reflection fraction R is defined as R = where fi is the solid 

angle subtended by the Compton-thick matter to the X-ray source; for a 



plane parallel slab, R = 1 . 



combined Chandra and VLBI observations imply for the kpc-scale 
jet bulk Lorentz factor a range between ~ 3.5 and 7.2. 

In the y-ray band, a 3 c r Egret detection was reported for this 
source JSguera et alj|2005t lHartman. Kadler. & T ueller 2008)), as 
well as an association with a Fermi source jAbdo et alJl201Ch . im- 
plying a broad-band SED reminiscent of a de-beamed blazar. Com- 
paring y-ray and multiwavelength prope rties for a sample of F ermi- 
detected and Fermi-undetected BLRGs. lKataoka et sug- 
gest that the GeV emission from 3C 1 1 1 is most likely dominated 
by the beamed radiation from the nuclear region of the relativistic 
jet. 



3 X-RAY OBSERVATIONS 
3.1 Suzaku 

Suzaku observed 3C 111 at the HXD (Hard X-ray Detector, 
Taka hashi e"t"aT1l2007t) nominal pointing position on 22 of August, 
2008 for a total exposure time of about 122 ks. The log of the X-ray 
observations is reported in Table [2] first part. We used the cleaned 
event files obtained from version 2 of the Suzaku pipeline process- 
ing. Standard screening criteria were used, namely, only events out- 
side the South Atlantic Anomaly (SAA) as well as with an Earth 
elevation angle (ELV) > 5° were retained, and Earth day-time el- 
evation angles (DYEJ3LV) > 20°. Furthermore, data within 256 s 
of the SAA passage were excluded and a cut-off rigidity of > 6 GV 

was adopted. 

XIS (X-ray Imaging Spectrometer. iKovama et ai1l2007l) data 
were selected in 3 x 3 and 5x5 editmodes using grades 0, 2, 3, 4, 6, 
and cleaned for hot and flickering pixels. The XIS spectra of 3C 111 
were extracted from circular regions centred on the sources of 2.9' 
radius. Background spectra were extracted from two regions, each 
with the same area of the main target region, offset from the main 
target and avoiding the calibration sources. The XIS response and 
ancillary response files were produced, using the latest calibration 
files available, with the ftools tasks xisrmfgen and xissimarfgen, re- 
spectively. 

We tested time- variability within the Suzaku observation, gen- 
erating source light curves in the soft and hard energy ranges with 
a binning time of 5760 sec, the orbital period of the satellite, to 
remove any possible modulation related to the orbital condition. 
Fig. [2 shows the soft and hard X-ray light curves obtained with 
the XISs, and their hardness ratios. The data from XIS0, XIS3, and 
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Figure 1. Background-subtracted light curves of 3C 111 from the Suzaku XIS0 (left panels), XIS3 (middle panels), and XIS 1 (right panels) observations, with 
a 5760 sec bin. Upper panels: 0.5 - 2keV; middle panels: 2 - lOkeV; lower panels: HR = [RATE2-10 — Rateo.5_2]/[rate2-io + RATEQ.5-2]. In each panel, we 
report the mean rate [HR] for the whole observation. 



XIS 1 are plotted on separate panels, from left to right . The source 
displays a hint of variability, with the x 1 probability of constancy 
less than 30% and 5% in the soft and hard band, respectively. How- 
ever, the HR light curves show a lack of significant spectral vari- 
ability. Therefore, in the following spectral analysis we considered 
the time-averaged data over the whole duration of the observation. 

After checking that the spectra were consistent, we coadded 
the two front illuminated (FI) CCD spectrfl along with the re- 
sponse files, to maximize the signal-to-noise ratio. The FI spec- 
trum was then fitted jointly with the back illuminated (BI, the XIS1) 
spectrum. 

The net FI [BI] source spectra were rebinned to 1024 channels 
and then grouped with a minimum of 200 [100] counts per bin. 
Data were included from 0.6 to 10 keV [8.5 keV], observed frame; 
the region around the Si K edge between 1.7 and 1.9 keV were 
discarded in both the FI and BI XIS, due to uncertainties in the 
calibration at this energy. The net background-subtracted source 
count rates for 3C 111 were 0.609 + 0.002, 0.712 ± 0.003 and 
0.676 ± 0.003 counts s~' for XIS0, XIS3 and XIS1 respectively, 
with a net exposure time after screening of 109 ks. 

For the HXD-PIN data reduction and analysis we followed 
the Suzaku data reduction guide (the ABC guidqj Version 2). We 
used the rev2 data, which include all four cluster units, and the best 
backgro und available, which account for non-X-ray background 
(NXb3 K okubun et ai1l2007l) . 

At the time of this analysis, two NXB files are available: back- 
ground A or "quick" background and background D or "tuned" 
background. We adopted the latter, which is the latest release and 
which suffers lower systematic uncertainties of about 1.3%, corre- 
sponding to about half uncertainty of the first release of the NXB. 
We confirmed this choice as most reliable estimate of the NXB by 
comparing the background A or D light curve to the light curve 
obtained from the Earth occulted data (from Earth elevation angles 
ELV < -5): indeed, the Earth occulted data give a representation 



XIS0 and XIS3, found to be consistent within 5% cross-normalization 
uncertainties 

3 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 

4 ftp://legacy.gsfc.nasa.gov/suzaku/doc/hxd/suzakumemo-2008-03.pdf 



of the actual NXB rate, as this neither includes a contribution from 
the source nor from the cosmic X-ray background. 

The source and background spectra were extracted within the 
common good time interval and the source spectrum was corrected 
for the detector deadtime. The contribution of the diffuse cosmic 
X-ray backg round counts was simulated using the spectral form of 
iBoldd « 1987), assuming the response matrix for diffuse emission, 
and then added to the NXB. With this choice of background, 
3C 111 is detected up to 70keV at a level of 19% above the 
background. The net exposure time after screening was 102 ks. 
The HXD-PIN spectrum was binned in order to have a signal- 
to-noise ratio greater then 10 in each bin, and the latest response 
file released by the instrumental team was used. The count rate 
in 14 - 70keV is 0.081 ± 0.002 counts s . Assuming a single 
cutoff power-law component (T ~ 1.6 and E c ~ 100 keV) this 
corresponds to a Fi^okeV ~ 4 x 10~"ergs cirr 2 s _1 . In the spectral 
analysis, we used a cross-calibration constant of 1.18 between 
the HXD and XIS spectra, as suggested by the Suzaku-HXD 
calibration team for observation at the HXD nominal pointing 
position^ 

Following the prescriptior0, we reprocessed the GSO data 
from the unscreened events, using the new gain calibration as of 
August 2010, in order to apply the correct GSO gain history file. 
Note that using the new GSO background instead of the one with 
the old gain calibration implies a loss of time of ~ 33%. As for the 
PIN, the GSO/NXB files was used in conjunction with the screened 
source events file to create a common good time interval. The cos- 
mic X-ray background, not included in the background event files, 
can be neglected, being less than 0.1% of the total background rate 
in the GSO. Both source and background spectra were corrected for 
the detector deadtime, and the former was rebinned for the back- 
ground subtraction, because the GSO background is created in 64 
bins. 

3C 111 is marginally detected with the HXD/GSO. The 
background- subtracted GSO count rate in the 50 - 200 keV band 
is 0.06 ± 0.03 counts s -1 , corresponding to a S/N ~ 1.3 (0.3% of 
the total), for a net exposure time of 68 ks. 



http://heasarc.nasa.gov/docs/suzaku/analysis/watchout.html 
http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/hxd_repro.html 
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Table 2. Observation log. 



Satellite 


Sequence N. 


Start (U.T.) 


Stop (U.T.) 


Detector 


Exposure 


Count Rate 












pes] 


[counts/s] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


Suzaku 


703034010 


2008-08-22 09:37:00 


2008-08-25 04:24:00 


XIS0 


108.9 


0.609 ± 0.002 










XIS3 


108.9 


0.712 ±0.003 










XIS1 


108.9 


0.676 ± 0.003 










PIN 


101.9 


0.081 ±0.002 










GSO 


68.1 


0.060 ± 0.023 


XMM-Newton 


0552180101 


2009-02-15 19:41:57 


2009-02-17 04:00:49 


pn 


59.1 


5.113 ±0.009 






2009-02-15 19:38:56 


2009-02-17 04:00:34 


MOS2 


60.9 


2.046 ± 0.006 






2009-02-15 18:59:02 


2009-02-17 04:10:44 


RGS1 


83.3 


0.086 ± 0.001 






2009-02-15 18:59:10 


2009-02-17 04:10:44 


RGS2 


83.2 


0.106 ±0.001 


XMM-Newton 


0065940101 


2001-03-14 12:56:20 


2001-03-15 01:20:36 


RGS1 


32.3 


0.132 ±0.002 






2001-03-14 12:56:20 


2001-03-15 01:20:36 


RGS2 


31.5 


0.162 ±0.002 



(1) X-ray satellite. (2) Sequence number of the observation. (3) Starting date (in year-month-day) and time (in hh:mm:ss) of the observation. (4) End 
date (in year-month-day) and time (in hh:mm:ss) of the observation. (5) Detector on-board the satellite. (6) Exposure time after screening was applied 
to the data. (7) Net source count rate after screening and background subtraction, as observed in the 0.6 - 10 keV (XIS0 and XIS3), 0.6 — 8.5 keV 
(XIS1), 14 - 70keV (PIN), 50 - 200 keV (GSO), 1 - lOkeV (pn, pattern 0), 0.4 - lOkeV (MOS2, pattern 0), and 0.65 - 2keV (RGS) energy ranges. 



3.2 XMM-Newton 

XMM-Newton observed 3C 111 twice, in March 2001 for about 
a total of 45ksec (Obs. ID 0065940101), and in February 2009 
for about a total of 120 ks (Obs. ID 05 52180101). The an alysis of 
the former observation is presented in [Lewis et al.l d2005h : the pn 
2.4 - lOkeV light curve from the se cond observation is inc luded 
in the X-ray monitoring presented bv lChatteriee et al". 

IHoH 

their 

"longlook data"). Here we present the analysis of the still unpub- 
lished X-ray spectra obtained in the second pointing. 

The observations were performed with the European Photon 
Imaging Camera (EPIC), the Optical Monitor (OM) and the Re- 
flection Grating Spectrometer (RGS). The observation details are 
reported in Table [2] second part. In this paper, we concentrate on 
the data in the X-ray band. 

The three EPIC cam eras (pn, MOS 1 , and MOS2; IStriider et all 
1200 ll ; iTurner et alj|200lh were operating in Small Window mode, 
with the thin filter applied. The XMM-Newton data have been pro- 
cessed, and event lists for the EPIC cameras were produced, using 
the Science Analysis Software (SAS version 10.0.2) with the most 
recent calibrations. 

EPIC event files have been filtered for high-background time 
intervals, following the standard method consisting in rejecting pe- 
riods of high count rate at energies > lOkeV. Given the X-ray 
brightness of the source, to avoid pile-up effects only events corre- 
sponding t o pattern have be en used (see the XMM-Newton Users' 
Handbook; iNess et aljlioiol) . We have also generated the spectral 
response matrices at the source position using the SAS tasks arfgen 
and rmfgen. 

Source counts were extracted from a circular region of radius 
25" and 30" for the MOS2 and pn cameras, respectively. For the 
pn camera, background counts were extracted from two source- 
free circular regions in the same chip of 50" radius each. Back- 
ground counts for MOS2 were extracted from similar regions but 
selected outside the central chip, symmetric with respect to the 
source position. Regarding the MOS1, in addition to two bad pixel 
columns near the source, a careful look to the images extracted 
in rawx,rawy coordinates revealed the presence of two more dark 
columns falling in the centre of the source, not visible using x,y 
coordinates, strongly affecting the spectral shape. Therefore, we 



decided to not use the data from the MOS1 camera, limiting the 
analysis of XMM-Newton only to the MOS2 and pn spectra. 

According to the result of the SAS task epatplot applied to 
the cleaned event file, pn data are affected by "X-ray loading", i.e. 
the inclusion of X-ray events in the offset map calculation^ This 
phenomenon, observed in very bright sources, is a consequence of 
excessive count rate, and occurs at lower source count rate than the 
pile-up threshold. The effects are a systematic energy shift to low 
energy and pattern migration from doubles to single events. As the 
pn offset maps are calculated on-board at the start of each expo- 
sure, for this camera the effects are hard to be reduced, therefore 
the pn data are more easily prone to this problem than the MOS 
one. At present, a reliable method to correct for this effect is not 
available. We then decided to discard the pn data in the soft energy 
range, more affected by the problem. The spectral shape of the pn 
spectrum above 1 keV is consistent with the MOS2 one, and with 
the Suzaku spectrum; this confirms that the selected band is not af- 
fected by the X-ray loading problem. Therefore, during the analysis 
we considered the energy range between 0.4 [1] and 10 keV for the 
MOS2 and the pn, respectively. 

As previously performed for Suzaku, we checked pn and 
MOS2 data for variability during the observation using the SAS 
task Icplot. No bin shows significant deviation from the mean value, 
allowing us to perform the spectral analysis over the time-averaged 
spectra. Both pn and MOS2 spectra were grouped to have at least 
200 counts in ea ch energy bin. 

The RGS (den Herd er et alj|200ll) data have been reduced us- 
ing the standard SAS task rgsproc, and the most recent calibration 
files; after filtering out the high-background time intervals, the total 
exposure times are ~ 83 ks for both RGS1 and RGS2. The RGS1 
and RGS2 spectra were binned at twice the instrument resolution, 
A/i = 0.2 A. 

3.3 Swift BAT 

The averaged BAT spectrum of this radiogalaxy was obtained from 
the 58-month survey archive (SWIFT J0418.3 + 3800). The data 
reduction and extraction procedure of the 8-channel spectrum is 

7 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0050-l-0.ps.gz 
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described in iBaumgartner et all j201dh . To fit the preprocessed, 
background- subtracted BAT spectra, we used the latest calibration 
response diagonal . rsp as of December 2010. 

3C 111 was detected in the 15-100 keV band with a count rate 
(1.02±0.03)x 10~ 3 counts s" 1 , which c orrespond to a 14—195 k eV 
flux of 1.2 ± 0.1 x 10-'° ergs cm" 2 s" 1 jBaumgartner et al.ll2010h . 

Fitting the averaged BAT spectrum with the HXD best-fit 
model (T ~ 1.6 and E c ~ 100 keV), we found F^ T 10keV ~ 7 x 
10" 11 ergs cnr 2 s -1 , while F fSnimv ~ ^ x " er 8 s cm 2 s~'( see 
Table O. The comparison between the BAT and the PIN spectra 
suggests a possible hard X-ray variability on longer timescales: 
fitting the two spectra together and allowing to vary the relative 
normalizations, we found that we need the cross-normalization be- 
tween the two instruments to be ~ 1.6, well above the uncertain- 
ties in the relative normalizations. The 58-month BAT lightcurv^l, 
provided with the catalogue, suggests the presence of hard X-ray 
variability and that the Suzaku observation took place during an 
upswing after reaching a minimum of the 14 - 195 keV flux. We 
caution that the Swift spectrum used here is the averaged spectrum. 



4 X-RAY SPECTRAL ANALYSIS 

All spectral fits to the X-ray data were performed using XSPEC 
v. 12.6.0q. The significance of adding free parameters to the model 
was evaluated with the F-test, with associated probability All 
uncertainties are quoted at the 90% confidence level for one pa- 
rameter of interest (A^ 2 = 2.71). Unless otherwise stated, the 
figures are in the rest-frame of 3C 111, and fit parameters are 
quoted in the same frame. All the models discussed in the follow- 
ing assume Galactic absorpti on with a column den sity of Ah. Gal ~ 
2.91 x 10 21 cnr 2 (see Tablefu lKaiberla et al.l2005h . It is worth not- 
ing that for 3C 111 the estimate of the total Galactic hydrogen col- 
umn density is not trivial: due to the presence of a dense molecular 
cloud along the line of sight, its value can as high as A'h, Gal ~ 
1.3 x 10 22 cnr 2 jLewis et al. 2005, and references therein), and it is 
expected to vary by several times 10 21 cnr 2 . Therefore we allowed 
the local absorption to vary up to 1.3 x 10 22 cnr 2 ; for this absorber 
we used the p habs model in XSPEC, adopting cross-sections and 
abundances of lWilms. Allen, & McCravl feOOof) . 

4.1 Suzaku 

The medium-energy spectrum of 3C 111 confirms the presence 
of strong iron K emission detected in previo u s obs e rvations 
.Lewis et alj 120051 ; iGrandi, Malaguti, & FiocchU 120061 : iDadinj 
l2007h . We first fitted the XIS and PIN data (between 0.6 and 
70keV) with a single absorbed power-law component, excluding 
the energy range where the iron K complex is expected (5-8 keV). 
After including the 5 - 8keV data, this model clearly provides a 
poor fit |^ 2 /d.o.f. = 1046.70/743; T = 1.61 + 0.01 and A H = 
(1.20 ± 0.02) x 10 22 cm" 2 ]; furthermore, large-amplitude residuals 
are present both in emission and in absorption (see Fig.|2j. No ad- 
ditional intrinsic A'h is required (Ay 2 /Ad.o.f. = 0.06/1); however, 
taking into account the strong variability of Ah, Gal for 3C 111, we 
cannot exclude that some of the observed absorption is at the source 
redshift. 

8 http://swift.gsfc.nasa.gov/docs/swift/results/bs58mon/mosaic_crab_lc/ 
BAT_58m_crabweighted_monthly_SWIFTJ0418.3+3800.1c.gz 

9 But see the caveats in using the F-test to m easure the significance of nar- 
row lines described in Protassov et al. (200-1). 
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Figure 2. Data-to-model ratio for the Suzaku spectra (red open squares, 
XIS1; black filled circles and line, FI XIS). The model is an absorbed 
power-law component (r ~ 1.61) fitted ignoring the 5 - 7.5 keV energy 
band, where the iron K complex is expected. For graphical purposes, the 
data have been binned to have a significance of 17cr. The positions of the 
iron K lines (Fe Ka at 6.4 keV and Fe K/? at 7.06 keV rest frame) are marked 
with a blue star. The vertical lines rep resent the centroids of the detected ab- 
sorption lines in the FI XIS data (see Tombesi et al. 20 10bh . 



The addition of a Gaussian line to the simple absorbed 
power-law component results in a large improvement in the fit 
(A^/Ad.o.f. = 244.5/3). The profile oftheFeKa line (E = 6.40± 
0.02 keV, a = 94 + _ 3 2 \ eV; F = 2.36+°^ x 10~ 5 photons cnr 2 s~\ 
EW = 94+jg eV) is marginally resolved, suggesting a possible con- 
tribution from the external disk and/or the BLR. However, con- 
straining the iron K line to be narrow (lOeV) does not leave strong 
residuals which could be due to an additional broad component. 
The strength of the Fe Ka line, coupled with the hard photon index 
(1.59 + 0.01) could suggest the presence of a reflection compo- 
nent. However, adding a refle ction component (parametrize d using 
the pexrav model in XSPEC; Magdziarz & Zdziarski 1995) to the 
absorbed cutoff power law plus iron line, we obtained a weak re- 
flection fraction, R = 0.16!°;^ Cf 2 /d.o.f. = 782.7/738). During the 
fit, the abundance was fixed to the Solar one; the photon indices of 
the cutoff power law and the reflection component, as well as their 
energy cutoff, were set to be equal. Allowing the high energy cutoff 
to vary we can set only a lower limit (E c > 75 keV). The inclinat ion 
angle (allowed to vary between 10° and 26°: lLewis et alj|2005l) . is 
also not constrained by the present data. We therefore decided to 
fix both parameters (see below). 

Although the energy centroid of the Fe Ka line suggests a 
low ionization reflector, we tested for a possible ionized reflector. 
We used an updated model for the Compton reflection off an op- 
tically thick photoionized slab of gas, whi ch includes the iron K 
emission line (reflionx m odel in XSPEC; IRoss. Fabian. & Yound 
ll999l : lRoss & Fabianll2005l) : as expected the ionization of the re- 
flector is found to be low, £ < 13.9erg cm s _1 , being the ionization 
parameter mainly driven by the energy centroid of the Fe Ka line. 

As the reflection component is weak and poorly constrained, 
we bounded its normalization to that of the emission line: to fit 
the two components in a consistent way, the ratio of their normal- 
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Table 3. X-ray fluxes and luminosities. 



3C 1 1 1 - Suzaku data - Baseline model 



Satellite 


[10 


Flux" 

-P -2 
ergs cm s 


- 1 ] 




0.5 - 2 keV 


2- lOkeV 


14-7()keV 


Suzaku 
XMM-Newton 
Swift 


9 cc+0.05 
-Q.D6 

6.83!™ 


48.32™ 


40 21 +IW 
71.19 ±2.45 


Satellite 




Luminosity'' 
[10 43 ergs s- [ ] 






0.5 - 2 keV 


2- lOkeV 


14-7()keV 


Suzaku 
XMM-Newton 
Swift 


J -^*-0.15 

1Z82$* 


11 03 +im 
11. uj „ 27 

27 54+0-62 


21 69 +1 

—1.03 

39.35 ± 1.36 



Errors are quoted at the 90% confidence level for 1 parameter of interest 
(AX 2 = 2.71). "Observed fluxes. * Intrinsic luminosities, corrected for 
local absorption, as well as for the effects of the photoionized gas. 



izations was assumed to be the expected onq 
illuminated b y a flat continuum (T ~ 1.6), 



| for a face-on slab 
e. N m ,,JNn 



1.2 x 1(T 2 ^George & Fabianl 1 199 it iNandra & Georgef 1 19941 ; 



iMatt. Fabian. & Rosd[l996l ; iMatt. Fabian. & Revnoldsll 19971) . The 
R parameter was fixed to 1, the inclination angle to 19° 
( lKadleretalJl2008l) . and the energy cutoff of the primary power 
law and of the reflection component to lOOkeV, a value consistent 
with the lower limit previously found. Higher values of E c imply a 
slight decrease in the quality of the fit. We also included a neutral 
Fe Kf3 line (at an energy fixed to 7.06 keV) with a flux o f 13.5% 
of the Fe Ka line dPalmeri et al.ll2003l ; lYaqoob et alj|2010h and the 
same width. 

The baseline model then consists of: a rather flat cutoff power 
law (r = 1.61 ± 0.01), and a reflection continuum+two narrow 
iron K lines. The width of the line, allowed to vary during the 
fit, is a = 70+jg eV. The strength of the reflection, estimated 
as the ratio of the pexrav and the direct power-law normaliza- 



tions, is iVpexrav/Ncutoffpi = 0.35 + 0.06. The column density is 
N H = (1.20 ± 0.02) x 10 22 cirr 2 . This model provides an ac- 
ceptable description for the broadband X-ray spectrum of 3C 111 
O^/d.o.f. = 793.8/740). 

However, some residuals suggest the presence of lines in emis- 
sion, around 0.9 keV, and in absorption, between 6 and 9 keV (see 
Fig. [3] and Fig. [4}. To account for this residuals we first added to 
the baseline model 3 Gaussian lines (2 absorption and one emis- 
sion lines). Including a narrow emission Gaussian component (a 
fixed to 10 eV) at low energies, the improvement in the fit is sig- 
nificant (A^ 2 /Ad.o.f. = 10.2/2, considering only 736 XIS bins be- 
tween 0.6 and lOkeV); the energy is E = 0.89 ± 0.01 keV, with 
EW = 32+ljjeV. Tentatively, the line can be identified with the 
O vm RRC or the Ne ix triplet lines. The two absorption lines 
are detected at the rest-frame energies of E = 7.23*^ keV (<r 
fixed to 50eV; EW = -15.2^° 9 8 eV) and E = 8.76 ± 0.20 keV 
(EW = -48.9! 2 ^ eV), with a A^ 2 /Ad.o.f. = 9.5/2 and 10.9/3, re- 
spectively. We refer to the Monte Carlo simulations presented in 
Tomb esi etai] £2010b) for a more robust estimate of the signifi- 
cance of the high-energy absorption lines in this Suzaku observa- 
tion of 3C 111. We note that in their analysis the authors considered 
only the FI data, and adopted for the 3.5 - 10.5 keV continuum a 
single power law, resulting in a flatter photon index T ~ 1 .47 and 
a higher continuum model above ~ 7 keV. However, given the lim- 



http://www.jca.umbc.edu/~george/html/science/seyferts/fgfplus.html 
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Figure 3. Data-to-model ratio for the Suzaku spectra below 1.5 keV; red 
open squares, XIS1; black filled circles and line, FI XIS. For graphical pur- 
poses, the data have been binned to have a significance of 10cr. The adopted 
model is the baseline model (an absorbed cutoff power law and a reflection 
component, with superimposed two narrow Gaussian lines, the Fe Ka at 
6.4 keV and the Fe Kfl at 7.06 keV rest frame). The centroid of the emission 
line when one more Gaussian is included in the model is marked with a blue 
star (emission line at E ~ 0.89 keV). 
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Figure 4. Data-to-model ratio for the Suzaku spectra above 5.5 keV; red 
open squares, XIS1; black filled circles and line, FI XIS. For graphical pur- 
poses, the data have been binned to have a significance of 17cr. The adopted 
model is the baseline model (an absorbed cutoff power law and a reflection 
component, with superimposed two narrow Gaussian lines, the Fe Ka at 
6.4 keV and the Fe KfS at 7.06 keV rest frame). The centroids of the absorp- 
tion lines when two more Gaussians are included in the model are marked 
with two vertical lines (absorption line at E ~ 7.2 keV and E ~ 8.8 keV). 



ited energy band considered, the effect of this continuum model 
does not strongly affect the detection significance of the absorption 
lines. If we consider only FI data, as the XIS1 is less reliable in 
the iron K energy range, the significance of the absorption lines is: 
A^ 2 /Ad.o.f. = 15.5/2 and 15.5/3 for the line at E ~ 7.23 keV and 
E ~ 8.76 keV, respectively. 

A detailed study of the absorption features and their interpre- 
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tation in terms of ultra-fast outflows r elated with the accretion disk 
is reported in iTombesi et al] teOlObh . Briefly, the most plausible 
association for the first absorption line is the resonant Fe xxvi Lya 
transition at 6.97 keV. If the line is identified with this atomic tran- 
sition, we have to assume a velocity shifts of v ou t 0.04 c. Con- 
cerning the second feature, it is likely a blend of different lines of 
the Fe xxvi Lyman series (namely, Ly/3, Ly-y, and Lyr5). In a more 
physically consistent modelling of the absorption feature, the ab- 
sorber is parametrized using a grid of photoio nized absorbers gen- 
erated by the XS TAlf"! photoionization code dKallman & Bautistal 
l200ll : lKallman et al]20 04). The photoionized gas is assumed to be 
fully covering the X-ray source with turbulence velocity v tllrb = 
500 km s , and illuminated by a power-law continuum from 
0.0136 to 13.6 keV with photon index T = 2. It is important to note 
that a single photoionized absorber can account for both lines. The 
results of the fit are: high ionization state, logf = 5.3*^ erg cm s" 1 , 
with ~ 6.97 x 10 23 cnr 2 (N'° n > 1.45 x 10 22 cnT 2 at the 90% 
level), and blueshift velocity v out = 0.038 ± 0.006 c (x 1 = 768.2 
for 735 degree of freedom). We note that the high ionization state 
implied by log<f ~ 5.3 erg cm s -1 does not introduce a signifi- 
cant curvature in the primary continuum below ~ 6keV. The out- 
flow velocity is somewhat similar to the source cosmological red- 
shift, z = 0.0485, and thus could suggests a co ntamination due 
to local or intermediate redshift absorbers (e.g., | McKernan et alj 
120031: iMcKernan. Yaqoob. & Revnoldsll2004l . l2005l) . However, this 
scenario is ruled out by the high column density i mplied for this 
absorbe r. We refer the reader to the discussion in ITombesi et al.l 
(2010b, sect. 5.1), where the authors investigate in details this pos- 
sibility and conclude that it is not feasible. The addition of the pho- 
toionized absorber improves the fit of Ay 2 /Ad.o.f. = 13.8/3. 

Observed fluxes and intrinsic luminosities are summarized in 
Table [3j the unfolded best fit to the data is shown in Fig.\5\(upper 
panel). The best-fit parameters are reported in the first part of Ta- 
ble [4] note that the inclusion of the warm absorber does not change 
significantly the continuum photon index or the emission lines pa- 
rameters. 

As a final test, taking into account the low reflection compo- 
nent and the flat continuum observed, we checked the presence of a 
second continuum emission due to the jet. To this end, we consid- 
ered also the GSO detection of 3C 111, that extends the source 
spectral energy distribution (SED) up to ~ 200 keV. Adding to 
the best fit model a second power law, with a flatter photon in- 
dex (Tje, = 1.51 ± 0.03), we obtained for the primary cutoff power 
law a photon index in better agreement with values expected from 
disk emission, T = 1.70+qJ^. We fixed the parameters of the lines, 
as well as the column density, the ionization state, and the velocity 
of the outflow, to the best-fit values previously found. This model 
provides a good description of the broadband X-ray spectrum of 
3C 111,^/d.o.f. = 846.6/768 (see Fig.|5] bottom panel). Despite 
the inclusion of the GSO data, the improvement in the fit due to 
the addition of this component is only A_^ 2 /Ad.o.f. = 3.1/2 (F-test 
probability Pp = 75.6%). We note that the jet dominates only above 
lOkeV, thus the possibility to accommodate a jet component in the 
best-fit model is strongly related to the inclusion of the GSO de- 
tection in the fit. Moreover, the relative normalizations of the two 
power laws are degenerate, making it difficult to spectrally distin- 
guish the true jet component. 



http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html 



3C 1 1 1 - Suzaku spectra 
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Figure 5. Suzaku XIS and HXD spectra of 3C 111, binned to have a signifi- 
cance of 3cr. In both panels, the black lines represent the total best-fit model; 
the underlying continuum is modelled with a cutoff power-law component 
(green dash-dot-dotted lines) plus a weak reflection component (magenta 
dash-dotted lines). The iron K line complex is composed by a Fe Ka and 
a Fe K/J features (blue dotted lines). In addition, an X-ray emission line 
at ~ 0.9 keV is also shown (orange dotted lines). In the bottom panel we 
show the model when a component due to the jet is considered (light blue 
line); also reported is the detection in the GSO, extending the spectrum up 
to200keV. 

4.2 XMM-Newton 

The 0.6 - lOkeV XMM-Newton emission shape of 3C 111 is 
in good agreement with the best fit for the broad band contin- 
uum of Suzaku, both in the broad-band slope (T ~ 1.60, N H ~ 
1.04 x 10 22 cirT 2 ) and in the detection of a strong Fe Ka emission 
line (see Fig.[6jl. The main difference is in the observed flux, a factor 
~ 2.5 higher than the Suzaku flux, see Table [3] This is not surpris- 
in g when the RXTE lo ng-term light curve is considered (see fig. 1 
in lChatteriee et al]|201 ll) : the Suzaku observation took place dur- 
ing a recovery of the emission after an historical low state, while 
the XMM-JVwton data corresponds to the average 2 - lOkeV flux 
level. 

Taking into account the Suzaku results and the evident pres- 
ence of iron K emission in the XMM-Newton spectra, we then 
investigated the presence of an emission line plus a reflection 
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Table 4. Best-fit parameters for the the Suzaku (XIS+PIN, 0.6 - 70 keV) and the XMM-Newton [MOS2+pn, 0.4(1) - 10 keV] spectra of 3C 111. 
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1 61+ 002 
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0.35 ± 0.06 
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0.89 ± 0.01 


74 ±27 
74 

lot 
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ioof 


19+ 005 
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7.06t 


89+53 
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The adopted model is a cutoff power law plus a reflection component, covered by a photoionized absorber, with the addition of the iron K complex; 
for the Suzaku data, a low-energy narrow emission line is also included. Local N H are (1.21 ± 0.02) X 10 22 cm" 2 and (1.04 ± 0.01) X 10 22 cm for the 
Suzaku and XMM-Newton spectra, respectively. " Parametrized with an XSTAR component obtained assuming a turbulence velocity v lur b = 500 km s~ 1 . 
b The ratio of the Fe Ka line and the reflection component normalizations was fixed to 0.012 (see text); the Fe K/? and Fe Ka ratio line was fixed to 
0.135; the Fe K/3 line width was coupled to the Fe Ko- line width. c For the XMM-Newton spectra, the highly ionized absorber is not required by the 
data; the inclusion of this component implies: N™ = 0.21^ ^ X 10 22 cnT 2 , logf = 2.80 + °" 22 erg cm s"\ and v out = 0.0010 ± 0.008 c. t Parameter 
fixed. 

(1) N\i of the ionized absorber. (2) Ionization parameter. (3) Outflow velocity. (4) Cutoff power-law and reflection component photon index. (5) Cutoff 
energy. (6) Reflection fraction (calculated as the ratio of the reflection component and the direct continuum). (7) Line energy. (8) Line width. (9) Line 
equivalent width. ( 1 0) x 1 and number of degree of freedom. 
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Figure 6. Data-to-model ratio for the XMM-Newton spectra (red open 
squares, MOS2; black filled circles and line, pn). The model is an absorbed 
power-law component (T ~ 1.60) fitted ignoring the 5 - 7.5 keV energy 
band, where the iron K complex is expected. The blue stars and the verti- 
cal lines mark the position of the features detected in the XIS data, as in 
Fig. |2] For demonstration purposes, the data have been binned to have a 
significance of 20cr. 



component, the latter parametrized again with the pexrav model 
in XSPEC, with abundance fixed to the Solar one, R parameter 
fixed to 1, inclination angle fixed to 19°, photon index set to 
be equal to the slope of the cutoff power law, and energy cut- 
off fixed to 100 keV. As expected due to the lack of data above 
lOkeV, the errors on the reflection component normalization are 
rather high. The Fe Ka line parameters are: E = 6.41 ± 0.03 keV, 



= 26+™ eV, 



F = 1.50™** X 10- 5 photons cm" 2 s- 1 , EW = 
24+] 2 eV (y 2 /d.o.f. = 1206.6/1143). As the reflection component 
is poorly constrained, we again bounded its normalization to that 



of the emission line. We also added a narrow Fe K/3, bounded 
to the Fe Ka as before. Thus we obtained a best-fit parameters 
consistent with the Suzaku one: cutoff power law and reflection 
continuum photon index T = 1.58 ± 0.01, reflection strength 



Nn 



pexrav IN cutoffpl 



19 +005 



, Nh = (1.04 ± 0.01) x 10 22 cm" 2 . The 



Fe Ka line parameters are: E = 6.41 ± 0.03 keV, cr = 89 + " eV, 



38^' eV with re- 



F = 2.34+jJ-^ x 10~ 5 photons cnT 2 s"\ EW 
spect to the observed continuum. In particular, we note a lack of an 
evident response of the line to the increase in the primary emission. 
Clearly, with only two observations we cannot determine the time 
lag, or compare with th e limit found with the RXTE monitoring by 
Chatt eriee et al.l d201 ll) ; as previously mentioned, the correlation 
found between X-ray flux and Fe Ka line intensity shows a lag less 
then 90 days, a result inconsistent with a torus origin for the bulk of 
the emission, although not with the BLR. We note that, despite the 
different and more complex model derived here, also in our case 
we found that the EW decreases with the increasing of the flux (see 
their fig. 8). 

This model provides an acceptable description for the broad- 
band X-ray spectrum of 3C 111 (r/d.o.f. = 1219.7/1144). Ob- 
served fluxes and intrinsic luminosities are summarized in Table[3j 
the unfolded best fit to the data is shown in Fig. [7] while the best-fit 
parameters are reported in Table[4] 

No clear residuals in absorption at ~ 7 - 8 keV are present. 
Adding an absorption line with energy and width fixed to the val- 
ues previously found, we obtained for the E = 7.23 keV and the 
E = 8.76 keV an upper limit to the EW of -8.5 eV and -1 1.1 eV, 
respectively, formally consistent within the errors with the Suzaku 
results. 

Finally, even though the previous model already provides a 
good description of the XMM-Newton data, taking into account 
the Suzaku results, we tested the XMM-Newton data for the pres- 
ence of possible outflowing gas. When the XSTAR table previously 
used is applied to the XMM-Newton data, the gas parameters best- 
fitting the Suzaku spectra result to be inadequate to describe the 
EPIC data 0f 2 /d.o.f. = 1234.1/1144): the lack of absorption fea- 
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Figure 7. XMM-Newton EPIC spectra of 3C 111. The black lines repre- 
sent the total best-fit model: the underlying continuum is modelled with 
a dominant cutoff power-law component (dash-dot-dot-dotted green line) 
plus a weak reflection component (dash-dotted magenta line). The iron K 
line complex is composed by a Fe Ka and a Fe KJ3 lines (blue dotted lines). 
For demonstration purposes, the data have been binned to have a signifi- 
cance of 20<t. 



ture at high energies, in spectra with S /N good enough to detect it, 
may suggests that the outflowing gas has varied. To further assess 
the possible variability of this highly ionized absorber, we fixed 
its ionization and velocity to the Suzaku best fit values. Allowing 
only the column density to vary does not improve the fit; the pa- 
rameter is unconstrained, with a 90% confidence interval between 
and N™ = 1.2 x 10 23 cm" 2 . This suggests again the variable na- 
ture of this disk wind. Allowing all the wind parameters to vary, a 

lower column density = 2.1^-g x 10 2 ' cnT 2 , lower ionization 

-l 



parameter logf = 2.80^25 erg cm s , and zero outflow velocity 
v om < 0.001c, are required Cr/d.o.f. = 1201.1/1141). The im- 
provement in the fit with respect to the model without absorber is 
Ax 2 = 18.6 for 3 degree of freedom less. This improvement is not 
due to single absorption features, which are unexpected for similar 
values of log £ and N™, but instead is due to a better parameterisa- 
tion of the spectral curvature due to the absorption. 



4.2.1 The RGS spectra 

In order to investigate the emission feature at 0.8 - 0.9 keV de- 
tected by Suzaku, we looked at the RGS data. To increase the 
statistics, we also reduced and analysed the RGS data obtained in 
March 2001; exposure times after cleaning and count rates are re- 
ported in Table [2] As anticipated, we rebinned the RGS spectra in 
constant wavelength bins at twice the instru ment spectr al resolu- 
tion (A/l = 0.2 A) and we used the C-statistic dCashll979h available 
in XSPEC for the spectral fit: with this choice of binning, there are 
fewer than 20 counts per resolution bin. 

A first inspection of the RGS1&2 data reveals that for 
both observations the spectra are background-dominated below 
~ 0.65 keV. We then considered in our analysis only data between 
0.65and2keV. 

The two observations show a rather similar spectral shape, al- 
though with a decrease in flux of a factor of ~ 1.5 (from the oldest 
to the most recent). In Fig.[8] upper and middle panels, we show the 
data/model ratio for the two observations when an absorbed power 
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Figure 8. RGS background-subtracted data-to-model ratio for an absorbed 
power law, in the energy range 0.68 - 1.05 keV. From top to bottom, Obs. 
0065940101, Obs. 0552180101, combined spectra. The blue star marks the 
position of the emission line at ~ 0.89 keV. 



law is applied, with T fixed to the value found for the EPIC data, 
1.58 (see Sect. [O}. 

We then combined the two observations (see Fig. [8] lower 
panel). Fitting the combined spectra with the same absorbed power 
law, we obtained a reasonable description of the RGS continuum 
(C/d.o.f. = 171.9/105 for 108 total bins). However, linelike resid- 
uals are present between 0.85 and 0.95 keV. 

We then added to the AGN baseline continuum an unre- 
solved line in emission, modelled with a Gaussian component; 
the improvement in the fit due to the addition of this component 
is AC/Ad.o.f. = 15.2/2. When allowing the width to vary, the 
line is still unresolved. The emission line is detected at an en- 
ergy of E = 0.890 ± 0.002 keV, with a flux of F = x 
10~ 4 photons cirT 2 S" 1 and an equivalent width of EW = 8 ± 1 eV. 
The line parameters from RGS1 and RGS2 are consistent within 
the errors. 

The most likely identification of this feature is the O vm RRC 
at 0.872 keV, although cannot be excluded a contribution due to 
the He-like triplet of Ne ix, at energie^ 0.905 keV, 0.9 1 5 keV, and 
0.922 keV. The quality of the data prevents us from further investi- 
gating the properties of the line. 



5 DISCUSSION 

In this paper we presented the analysis of our Suzaku observation 
of 3C 111, as well as archival unpublished XMM-Newton spectra 
of the BLRG. 



http://physics.nist.gov 
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The Suzaku and XMM-Newton emission of 3C 111 is char- 
acterized by a hard continuum, showing weak reprocessing fea- 
tures (iron K line and reflection component). The continuum is 
rather flat, T ~ 1. 6, compared to the values typically found for 
RQ (T ~ 1.9; e.g. iReeves & Turned l200d : IPiconcelli et"aT]|2005l; 
Mateos et alj|2010h and BLRGs (T ~ 1.7 ; e.g. lZdziarski & Grandil 



200 ll ; iGrandi, Malaguti. & Fiocchill2006l) . The broad-band shape 
observed by Suzaku and XMM-Newton is in good agreement; the 
0.5-10 keV flux changes of a factor 2.5 in the 6 months separating 
the two observations. Such a strong flux variability is not surpris- 
ing in this source. The power spectral density anal ysis of the RXTE 
monito ring between 2.4 and lOkeV presented by IChatteriee et all 
d201 ll) reveals a break timescale of ~ 13 days; flares as strong as a 
factor of 2 - 3 in flux are observed in ~ 1 month. Previous XMM- 
Newton and RXTE observations caught the source in a higher flux 
state than Suzaku (by a factor of 3 in the 2 - lOkeV energy range) 
as well as in the new XMM-Newton data. Again, the X-ray spec- 
tral shape does no t change significantly: T ~ 1.63 - 1.75, R < 0.3 
dLewis et alj|2005h . Variations in the intensity of the emission are 
observed up to at least ~ 100 keV, as implied by the comparison 
between the HXD-PIN emission and the BAT averaged spectrum. 
Despite the continuum variation, the fluxes of the 6.4 keV Fe Ka 
line, clearly detected in both XMM-Newton and Suzaku datasets, 
are consistent. 



Table 5. Spectral indices and fluxes from the combined fit to the X-ray 
SED. 



Satellite (Year) 


r 


Flux" [10- 11 


ergs cm 2 s 1 ] 






0.5 - lOkeV 


10-200keV 


BeppoSAX (1998) 


1.67 ±0.06 


2.84 


7.42 


XMM-Newton (2001) 


1.68 ±0.01 


7.43 


15.78 


Swift (58-month) 


1.57 ±0.09 


3.90 


11.41 


Suzaku (2008) 


1.61 ±0.01 


2.22 


6.89 


INTEGRAL (2008) 


95+ - 43 


1.08 


12.35 


XMM-Newton (2009) 


1.61 ±0.01 


5.45 


14.05 



Errors are quoted at the 90% confidence level for 1 parameter of interest 



(Ax = 2.71). a Observed fluxes. 



0.002 keV. The measured fluxes are 

+0.71 v, i n-4 „w„„„ ™-2 „-l 



sistent energy, E = 0.890 

consistent within the errors: 1.43^2'Z] x 10 -4 photons cm 1 s~' and 
2.09+55 85 x 10 ~ 4 photons cm" 2 s" 1 for XMM-Newton and Suzaku, 
respectively. Similar features, produced by photoionized gas, have 
bee n detected in two other FRII radio galaxies, the BLRG 3C 445 
(e.g. lGrandi et aL 2007; Sa mbruna, Reeves. & Braitol2007l) and the 
NLRG 3C 234 jPiconcelli et alj|2008l) . " 



5.1 Emission and absorption features 

In addition to the Fe Ka-Fe K/3 complex, the most significant fea- 
tures detected in the XIS spectra are two absorption lines at E ~ 
7.24 and ~ 8.77 keV, respectively. The presence of such absorp- 
tion features are not unique to 3C 111: an analysis of 5 nearby, 
X-ray bright BLRGs s hows that these absorp tion lines are detected 
in 60% of the sources jTombesi et alj20 10b. where a detailed anal- 
ysis of these features is discussed). Their likely interpretation as 
blue-shifted iron K lines implies an origin from highly ionized gas 
outflowing with v out ~ 0.04 c, probably related with accretion disk 
winds/outflows. From the ionization parameter £ = L im /(nR 2 ) and 
the column density of the outflowing gas, under the reasonable as- 
sumption of thickness of the clouds AR = Nu/n lower than the 
distance, AR <K R, we can estimate the launch radius. Adopting the 
absorption-corrected ionizing luminosity between 0.0136 keV and 
13.6 keV as obtained from the Suzaku data (2.5 x 10 44 ergs s -1 ), we 
found a sub-parsec distance from the central BH, R < 9 x 10 15 cm. 

Assuming the m ean black hole mass estimated by 
IChatteriee et al.1 < l201ll) . M BH = 1.8 x 10 8 M o , the escape 
velocity at this radius is a factor 2 higher than the outflow velocity, 
v csc ~ 0.077 c. The hig hest black hole mass estimated for 3C 111, 
M BH = 3.6 x 10 9 M G iMarchesini. Celotti. & Ferraresg|2004 see 
Table [TJ, would imply an even worse situation, with an escape 
velocity a factor 10 higher than v out . The outflow is then unable to 
escape from the central region where is produced, and the absorber 
is most likely unstable. Indeed, at the time of the XMM-Newton 
observation the outflowing gas has varied. No absorption features 
are detected in the EPIC spectra. Testing a photoionized absorber 
as seen in the Suzaku data, we found a steady absorber with lower 
column density and lower ionization parameter. 

Finally, the XIS spectra also show an emission feature at ~ 
0.9 keV; the most commonly observed features at energies around 
E = 0.89 + 0.01 keV are the O vra RRC at 0.872 keV or the 
Neix triplet lines at 0.905 keV, 0.915 keV, and 0.922 keV; at the 
XIS resolution these features are unresolved. An emission line is 
also detected in the combined RGS spectra of 3C 111, with con- 



5.2 From X-ray to y-ray 

The comparison presented here between HXD-PIN emission and 
the BAT averaged spectrum at high energies, and XMM-Newton 
and XIS data at lower energies, confirms the strong variability of 
3C 111, which could be tentatively ascribed to the jet variability. 
Only broad-band monitoring observations will allow us to con- 
firm or rule out this possibility. No clear roll-over is detected in 
the Suzaku data, which can be explained by the presence of such 
a jet component. Indeed, if the GSO detection is robust, there is 
a marginal evidence of the emerging of a jet component at ener- 
gies higher than lOkeV, although the present data do not allow us 
to firmly establish its relative contribution with respect to the pri- 
mary continuum. The strong variability observed in the X-ray flux 
together with the flat photon index of the continuum, could suggest 
dilution of the disk emission due to a variable jet contribution. In 
principle, data with higher spatial resolution could allow us to in- 
vestigate the jet contribution below lOkeV, resolving some knots 
that are blended for XMM-Newton (and even more for Suzaku); in- 
deed, X-ray jet emission may arise on a range of scales, as the one- 
sided radio jet is detected at mas to arc-minutes. However, images 
from a Chandra observation l lHogan et alj|201 lh show that all but 
one the visible knots along the jet fall outside the extraction region 
of the XMM-Newton spectra, and none of them is bright enough 
to be detected in the automatic source detection and analysis per- 
formed with XASSISlF^l The lack of unresolved bright knots in the 
XMM-Newton spectra, that correspond to a higher flux level, make 
unlikely that the variation is jet-related. 

In the y-ray band, an Egret identification of 3C 111 as a 
y-ray source has been suggested, although with some debate on 
the real contributi on of 3C 111 to the total y-ray emission of the 
Egret counterpart dSguera et al.ll2005l ; lHartman. Kadler. & T ueller 
2008). Re cently, an as s ociatio n with a Fermi source has been re- 
ported by lAbdo et all d2010h . The fluxes detected in the y-ray 
band are: Fo.75-3McV < 4.2 x 10~ 10 ergs cirT 2 s and F3-30MeV < 

13 http://xassist.pha.jhu.edu/zope/xassist 
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3C 1 1 1 - Galactic Abs. -corrected 




Figure 9. Rest-frame high-energy SED of 3C 111. Red symbols are Suzaku spectra (XIS03, PIN, and GSO); above 10 18 Hz, green and purple circles are BAT 
and INTEGRAL data, respectively. BeppoSAX MECS and PDS spectra are plotted in light blue; orange symbols are XMM-N ewton MPS data from 2001, 
while blue spectrum is the XMM-Newton MOS2 spectrum analysed in this paper. Black filled squares are taken from literature (Hartman, Kadl er. & TuelleJ 
120081 : lAbdo et alj riblQj): from low to high energies, ASCA, CGRO (OSSE, COMPTEL, and Egret instruments), and Fermi data. Data are corrected for local 
absorption. Dashed black lines are the components of the best-fit model for the Suzaku data when the jet component is included; the solid line corresponds to 
the total model, with a cut off at ~ 10 MeV applied. 



4.3 x 10 10 ergs cm 
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[Egret, b oth on-board the Compton Gamm a 



Ray Observatory (CGRO); lHartman. Kadler. & Tuellerl 
fi-ioooev ~ 6x 10" 12 ergs cm" 2 s" 1 (Fermi: lAbdo et al] 



2008]; 



2010) 



These fluxes are highly above the value predicted by the best-fit 
model to the 0.6 - 70 keV Suzaku data without a jet component. 
This supports the hypot hesis of high-e nergy emission due to the 
jet in 3C 111. Recently, Katao kaet al.l d201lh ascribed the Fermi 
detection to a pc-scale relativistic emission from the jet. Extrapo- 
lating the dual power law model, including a possible jet compo- 
nent (see Sect. 14.1b , to high energies, we found for the jet com- 



ponent f J 3 ct 30MeV = 4.9 x 10" 10 ergs cm" 2 s" 1 . Taking into account 
the typical variability of jet emission, this flux is in substantial 
agreement with the upper limits obtained in the COMPTEL energy 
range. Qualitatively, a peak at ~ 100 MeV would cause the ob- 
served jet component to be compatible with the observed emission 
above 0. 1 GeV. 

This is evident looking at Fig. [9] where we present the high- 
energy historical SED for 3C 111 from X-ray to y-ray energies. 
We used the Suzaku, XMM-Newton, and Swift observations anal- 
ysed in this paper together with ASCA, XMM-Newton MOS (from 
2001), BeppoSAX (from 1998), INTEGRAL (from 2008), CGRO, 



and Fermi data. For safe of clarity, in Fig.[9]we do not show XIS1 
and pn (from 2009) data. We note that the Swift spectrum is av- 
eraged over more than 4 years. The strong variability previously 
mentioned is clearly evident when all the datasets are plotted to- 
gether. In order to understand the driver of the variability, we fitted 
all the datasets together, adopting the best-fit continuum model of 
the Suzaku data, i.e. an absorbed cutoff power law plus a reflection 
component and the iron K complex; since the comparison between 
Suzaku and XMM-Newton spectra shows that the lines do not vary 
with the continuum, we fixed the reprocessing features to the best 
fit value found for Suzaku. 

As a first step, we allowed only the relative normalization of 
the primary power law to vary between the observations. For ob- 
servations extending down to a few keV, also the column density 
was allowed to vary. This fit already provides a reasonable repre- 
sentation of the SED (x 2 /d.o.f. = 3177.6/2429), our aim being to 
investigate the variations in the main spectral properties (continuum 
slope and emission intensity), and not to find a detailed description 
of all the spectra. However, to further investigate the observed vari- 
ability, we also allowed the power-law photon index to vary; the ob- 
served change of intensity and spectral shape is summarized in Ta- 
ble[5]0f 2 /d.o.f. = 2972.8/2423). We obtained rather similar values 
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x 10-" ergs cm- 2 s" 1 (our XIS 
7.4 x 10"" ergs crrT 2 s" 1 (XMM-Newton from 2001). 



of the photon index, spanning a range from T = 1.57 ± 0.09 (Swift) 
to T = 1.68 ± 0.01 (XMM-Newton from 2001). The only outlier is 
the INTEGRAL observation, for which we obtain T = 0.95^. 
This is probably due to the low quality of the data, confirmed 
by the poor constraints to the value of T, coupled with the high- 
energy cutoff applied, E c = 100 keV. Fitting the data with a simple 
power law would increase a little the best-fit value of the photon 
index, although with strong residuals at the high-energy end of the 
spectrum. The column density varies between ~ 1.08 x 10 22 cm~ 2 
(XMM-Newton from 2009) and ~ 1.50 x 10 22 crrT 2 (BeppoSAX 
from 1998). As noted before, we observe significant variability in 
the 0.5 - lOkeV flux, from ~ 2 
data) to 

At higher energies, the stronger variation is observed between the 
Suzaku and the XMM-Newton (from 2001) data. 

This analysis shows that a brightening of the source does not 
imply a flatter power law, thus suggesting that the main driver of 
the increase in the flux may not be a strong increase in the jet com- 
ponent. Although our analysis is limited by the number of spectra, 
we note that the photon indices are well constrained thanks to the 
use of soft X-ray data as well as data above lOkeV. After exclud- 
ing the INTEGRAL point, the change in slope with the flux found 
here, adopting a more complex model, is marginally consistent w ith 
that obtained from the RXTE monitoring (Chatte riee et al.l201 ll see 
their fig. 8), which shows a similar trend albeit with a larger Ar and 
covering a large range in the 2 - lOkeV flux. 

In our observational program, we have BLRGs showing a 
spread of the jet angle with respect to the line of sight. Despite the 
strong uncertainties in the orientation, the closest candidate to be 
comp ared to 3C 1 1 1 i s 3C 3 82. Recently, from high-quality Suzaku 
data lSambruna et all d201ll) estimated for this BLRG a disk incli- 
nation of ~ 25° - 30°; at bigger scales, radio observations found a 
lower limit to the inclination angle of ~ 15° dEracleous & H alpern 
1998). The best estimate of the jet inclination for 3C 111 is ~ 19° , 
with a proposed range of 10° < ( < 26° ( Kad ler et a i] |2008l : 
iLewis et"al]|2005l) . Despite this similarity, the X-ray emission of 
the two sources is remarkably different, with 3C 382 showing an 
X-ray spectrum more similar to that observed from Seyfert galax- 
ies. This implies that the relative importance of the jet is not related 
only to geometrical properties (i.e., orientation with respect to the 
line of sight). Different velocities, or a bent jet, can play an impor- 
tant role. 

Finally, the X-ray emission of 3C 111 appears significantly 
different from that of a typical Seyfert galaxy in many aspects, e.g. 
flatter intrinsic emission and weak reflection features. On the other 
hand, its time variability shows properties that correlate with the 
black hole mass, following the same scaling law observed in stellar- 
mass black hole X-ray binaries and Seyfert galaxies, suggesting 
that a similar accretion pro cess is powering these different sys- 
tems dChatteriee et al.l201 lh . It is worth noting the large dispersion 
of int rinsic X-ray properties shown by RL AGN dSambruna et al.l 
120091 and references therein), similar to the one found out for the 
RQ sources. The overlapping of the distributions of photon in- 
dex and reflection features observed in both classes supports the 
idea of common accretion structure, with a second parameter in 
addition to the accretion rate that determines the jet production 
efficiency. The most likely candidate is the black hole spin, al- 
though how does it works in producing the observed RLs-RQs 
dichotomy, or the FRIs-FRIIs division, is still strongly debated 
(e.g., the "spin paradigm" coup led with intermittent je t activity, 
Sikora. Stawarz. & Lasotall2007l; or the "ga p paradigm", Garofalo 



20091 ; iGarofalo. Evans. & Sambrunal [2010). In particular, in the 



heuristic scenario proposed by Garofalo et al., collimated jets cou- 
pled with intermediate accretion efficiency, as observed in 3C 111 
(see Table [TJ, can occur for retrograde systems. This configura- 
tion also implies a larger size of the gap region between the in- 
ner edge of accretion disks and the black hole horizon, size that is 
connected to the energetics of the disk itself. The weakness of re- 
flection features is a natural consequence of such a configuration; 
moreover, larger gap regions can limit the presence of disk winds 
dKuncic & Bickn ell 2004, 2007). Unstable outflows, weaker than 
in RQ AGN, could be theref ore compatible with H igh Excitation 
FRII states, such in 3C 111 dButtiglione et alj|2010l) . while at the 
extreme end, radiatively inefficient accretion flows in Low Excita- 
tion Radio Galaxies (L/Leu ~ 10~ <5 ~ 7) ) might inhibit the launch of 
disk winds. 



6 SUMMARY 

In this paper we presented the X-ray spectra of 3C 111, obtained 
with Suzaku as part of our observational program devoted to study 
the high-energy emission of the brighter BLRGs. A recent XMM- 
Newton observation of ~ 120ksec is analysed as well. 

The 3C 111 emission extends up to ~ 200 keV, with a rather 
flat continuum; in spite of this, no signatures of jet emission are vis- 
ible below lOkeV. Suzaku observed the source at a minimum flux 
level, as shown by the RXTE monitoring. In the 6 months separat- 
ing the Suzaku and XMM-Newton observations, the flux increases 
by a factor of 2.5, a level of variability which is not unexpected in 
this source. We confirm the weakness of the reflection features, as 
found in previous observations. An iron K complex is clearly de- 
tected in both datasets, with a rather low EW. The intensities of the 
line detected by XMM-Newton and Suzaku are consistent, with a 
lack of immediate response of the line properties to the continuum 
variation. In addition to the iron K complex, absorption features due 
to a highly ionized, ultra-fast, nuclear outflowing gas are detected 
in the XIS data. At the time of the XMM-JVewton observation, the 
gas has varied: adopting the photoionized absorber as seen in the 
Suzaku data, then a lower column density and lower ionization state 
are required. Both XIS data and RGS spectra show the presence of 
a narrow line at ~ 0.89 keV, which is most likely identified with 
either the O vm RRC or the Ne ix triplet lines. 

The lack of detection of a roll-over in the primary emission 
is probably due to the appearance of the jet as a dominant compo- 
nent in the hard X-ray band, as suggested by the detection with the 
GSO on-board Suzaku above ~ 100 keV. From a qualitative point 
of view, the emission observed by Egret and Fermi above 0. 1 GeV 
is compatible with such a jet component assuming that this peaks 
at energies of ~ 100 Me V. If the detection is confirmed, the emer- 
gence of a jet component can explain differences observed in the 
high-energy emission of BLRGs similar in other aspects. The rel- 
ative importance of the jet component is not simply related to the 
system inclination, but can be associated to the velocity of the jet; 
a change of the opening angle of the jet passing from the nuclear 
region, where the observed X-ray emission is produced, to the most 
distant regions, where the radio emission is observed, may also play 
a role. 

Finally, we presented the high-energy historical SED for 
3C 111, from X-ray to y-ray. Our qualitative analysis suggests that 
the strong variability observed for 3C 111 is probably driven by a 
change in the primary continuum. This implies that simultaneous 
broad-band X-ray and y-ray monitoring is needed to unambigu- 
ously constrain the parameters of the disk-jet system, and correctly 
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understand the contribution of the different components (i.e., disk 
emission, reflection, jet) to the total emission. 
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